Interactions and screening in gated bilayer graphene nanoribbons 
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The effects of Coulomb interactions on the electronic properties of bilayer graphene nanoribbons 
(BGNs) covered by a gate electrode are studied theoretically. The electron density distribution and 
the potential profile are calculated self-consistently within the Hartree approximation. A comparison 
to their single-particle counterparts reveals the effects of interactions and screening. Due to the finite 
width of the nanoribbon in combination with electronic repulsion, the gate-induced electrons tend to 
accumulate along the BGN edges where the potential assumes a sharp triangular shape. This has a 
profound effect on the energy gap between electron and hole bands, which depends nonmonotonously 
on the gate voltage and collapses at intermediate electric fields. We interpret this behavior in terms 
of interaction-induced warping of the energy dispersion. 

PACS numbers: 81.05.Uw, 73.23.-b, 73.21.Ac 



I. INTRODUCTION 



Graphene and its bilayers are extremely good con- 
ductors with mobilities up to 20m 2 /Vs at room 
temperature^! 2 - which makes them interesting as ac- 
tive channels in field-effect transistors (FETs). These 
materials are also considered as building blocks for 
quantum electronics such as quantum computation and 
pscudospintronic devices To facilitate these applica- 
tions, large on-off resistance ratios are necessary^ Sev- 
eral mechanisms for opening up an energy gap in two- 
dimensional (2D) monolayer graphene have been pro- 
posed to achieve this goali 7 - It has been suggested that the 
application of an electric field between the two graphene 
layers breaks the inversion symmetry and induces an 
energy gap between electron and hole bands.— Oost- 
inga et al<2 have recently been able to detect such an 
electric-field induced energy gap in a graphene bilayer 
system. Infrared absorption spectroscopy furthermore 
demonstrated a gate-controlled band gap with gap en- 
ergies up to 250 mcV, which is in accordance with the 
self-consistent tight-binding calculations ! 10 ' 11 The size of 
the gap is determined by the electric field between the 
two layers which can be tuned in single— or double- 
gate geometries. Recently, the influence of a potential 
applied to a single gate on bilayer graphene has been 
studied theoretically within the continuum model, reveal- 
ing a roughly linear dependence of the gap width as the 
carrier density increases^ However, a nonmonotonous 
evolution of the gap width with the electron density has 
been predicted for graphene triple and quadruple layer 
systems, which has been attributed to trigonal warping of 
the band structure! 15 ' 16 Moreover, ab-initio density func- 
tional theory calculations have been used extensively to 
investigate the electronic structure of bilayer graphene 
and essentially confirm the behavior suggested by the 
tight-binding and continuum models*^ - — 



Due to the presence of edges, graphene nanoribbons 
(GNRs) reveal a much richer phenomenology than 2D 
graphene sheets 3"— GNRs possess energy gaps due to 
size quantization and edge states are formed with prop- 
erties depending on the width and the type of edgesi 23 ' 24 
Edge disorder may result in transport gaps around the 
neutrality point £!L2ki2L For gated GNRs, the electronic 
structure of monolayer systems has been modeled by 
Fernandcz-Rossier with the image charge method in 
the Hartree approximation^ The inhomogeneous charge 
density and potential across the nanoribbons were found. 
Also, the small transverse size has a strong impact on 
the classical and quantum capacitances of monolayer 
GNRs.— 



It is thus self-evident that interactions and imperfect 
screening should influence the electronic properties in bi- 
layer GNRs (which we will refer to as BGNs below) as 
well and may modify the formation of the gate- voltage in- 
duced energy gap. However, neither the gate electrostat- 
ics, nor the self-consistent band structure and the band 
gap formation in BGNs (in contrast to bulk bilayers) have 
been studied before. In the present paper, we address 
this issue by taking Coulomb interactions in BGNs into 
account and studying the electron density distribution, 
the Hartree potential and the energy dispersion, focus- 
ing particularly on the evolution of the electric field in- 
duced band gap which can be regarded as a single- valued 
parameter for the relevance of the interaction effects. 



The paper is organized as follows. In Sec. II we intro- 
duce the structure to be investigated and formulate our 
model. The results are presented and discussed in Sec. 
III. Sec. IV contains a summary and conclusions. 
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II. STRUCTURE AND MODEL that the type of edge of the BGN has only a marginal 

influence on our results and that finite-lengths effects are 
Throughout this work we consider long BGNs stacked neglected. The tight-binding Hamiltonian reads^ 
in the Bernal form with ideal armchair edges. We note 



j) -7o ( a Zfr>i + h - c - 



7i ^fat.^,. + h.c.) - 73 Y ( 5 M fl 2.i + h - c -) 

1 



(1) 



r 



(a) 



gate 



Z 



diel ;ct ic 



+0 

P: 



t8 BGN pT 




FIG. 1: (a) Scheme of the BGN device. The arrows indicate 
the direction of electric field for V g > 0.(b) Illustration of the 
band structure shift induced by the gate voltage V g , and of 
the voltage components V c and V q . 



where (b^) are the creation operators at sublatticc 
A (B) in layer I = 1,2 at site r*. We adopt the com- 
mon graphite nomenclature and include the three cou- 
pling energies 70, 71 and 73, where 70 = 3.16 eV rep- 
resents the intra-layer nearest-neighbor coupling energy, 
71 = 0.39 eV is the coupling energy between sublatticc 
B and A' in different graphene layers, and 73 = 0.315 cV 
the hopping energy between sublattice A and B' in the 
lower and upper layers, respectively. Vi is the the on-site 
potential arising from the Coulomb interaction of the in- 
duced charges the density of which is to be calculated. 

In the capacitor structure shown in Fig. QJa), a BGN 
of width w in y-direction and layer spacing 6 = 0.335 nm 
in z-dircction extends along the x-direction and is em- 
bedded in an insulating medium (we assume a relative 
permittivity e r =3.9 adequate for Si02). The structure 
is covered by a homogeneous metallic gate extending in 
the xy - plane at distance d. Electrical access is enabled 
by connecting a voltage source to the metal gate and the 
lower BGN layer. A neutral BGN is obtained by setting 
the on-site potential to zero, corresponding to the Fermi 
energy Ep = 0. As a gate voltage V g is applied, excess 



charges pi = etii are induced in the lower (i=l) and up- 
per (i=2) graphene layer summing up to the total charge 
density p = en. Here, n, denotes the carrier density in 
layer i and e the electron charge. This is balanced by 
the charge density a of opposite polarity at the metal- 
oxide interface such that the Fermi energy departs from 
the neutrality point, while the charge neutrality of the 
whole system is conserved. In the following calculations, 
we subtract V g from both the gate electrode as well as 
from the BGN which does not modify the electrostatics 
of the problem. This way, the gate metal-insulator inter- 
face is kept grounded which allows us to take advantage 
of the image charge technique creating a charge density 
at a distance d above the metal-insulator interfaced 

Within our framework, eV g is defined as the difference 
between the electrochemical potential u and the neutral- 
ity point of the discharged BGN u , 



eV g = u - u (h 



(2) 



see Fig. [T] (b) for an illustration. The displacement of 
the chemical potential is composed of the classical term 
eV c associated with the electrostatic potential, and of the 
quantum term eV q describing the accumulation of extra 
charges in the energy bands, 



v g = v c 



v q - 



(3) 



The charge distribution screens the lateral electric field 
components at the gate. Taking the Coulomb interac- 
tions between the charges and their images into account, 
the Hartrce potential at the site r in the continuous limit 
can be computed as 



I 

with ^ = f and for the lower and the upper layer of the tion between the charges in each layer and their corre- 
BGN, respectively. The first term describes the interac- 
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sponding mirror charges, while the second term accounts 
for the interaction of the charges in one layer with the 
mirror charges related to the other layer. 

In order to determine the charge density, the local den- 
sity of states D(r, E) (LDOS) at site r and energy E is 
calculated with the help of the real-space Green's func- 
tions technique described in Ref. [35|. The electron den- 
sity at site r is then given by 



n(r) 



D sc (E)f{E)dE 



D (E)f(E)dE (5) 



where f(E) is the Fermi-Dirac distribution. D SC (E) and 
Dq(E) represent the LDOS at energy E with and without 
including the self-consistent interactions, respectively In 
the actual calculations, the second term is replaced by 
4/ v / 3a 2 with a = 0.246 nm which represents the positive 
charge background of the ions. Since the Hartrce poten- 
tial, Eq. (@| depends on n(r) which is a solution of the 
Schrodinger equation with the Hamiltonian given by Eq. 



([T]), these equations are solved by the Broyden iterative 
method. 

Some insight into the basic physics of gated BGNs is 
obtained from the analytical expression for the potential 
difference between the layers induced by gate voltage. 
Since 5 <C d, we can approximate the BGN by a con- 
ducting strip with the total charge density p. The rela- 
tion between gate voltage and the charge density for the 
system formed between the conducting strip of width w 
located at the distance d apart from a semi-infinite plate 
is given by 3 -^ 



TT£oE r 



2d arctan 




The potential difference between the graphenc layers can 
be calculated by integrating the electric field in between, 
i.e., AV = fj +S E(z)dz. By using the method of image 
charges, one finds that E(z) is determined by 



E{z) = Ef^iz + d) + E%*{z + (d + 6))+ E^(z - d) + E^(z -(d + S)) 



(7) 



r 



where the first two terms in the summation correspond 
to the image charges. The intensity of the electric field 
at distance z from the middle of a strip is given by 3 -^ 
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Performing the integration and using Eqs. (JS])-©, we 

get 



S arctan {^)V g 



2d arctan i 
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(9) 



where Ap = p-2 — p\ is the charge density difference. In 
the derivation of Eq. (|9]), a homogeneous distribution of 
the charge density in the BGN has been assumed. This 
assumption is a reasonable approximation for low charge 
densities. 



III. RESULTS AND DISCUSSION 

Figs. [Ha) and (b) show the electron density distribu- 
tions and Hartree potentials of BGN with N = 49 sites 
in transverse direction, corresponding to w = 12nm for 
different gate voltages. The calculations are performed 
over one unit cell containing 47V sites. To clarify the 
role of electron-electron interactions we show the self- 
consistent Hartrce (right parts) and non-interacting (left 
parts) calculations. In our context, non-interacting cal- 
culations correspond to the Hamiltonian ([1} where the 



potential difference between the layers is given by Eq. 
(|9]) with Ap set to 0. Both non-interacting and self- 
consistent calculations show pronounced oscillations be- 
tween neighboring sites. Because of this we also plot for 
clarity the density distributions averaged over six neigh- 
boring sites. Within the self-consistent treatment, the 
electrons tend to accumulate along the edges as a man- 
ifestation of electronic repulsion. As the gate voltage is 
increased, the edge accumulation becomes more promi- 
nent. It is worth pointing out that it is primarily the 
applied gate voltage which determines this distribution, 
resembling that one found in monolayer GNRsi 30 ' 31 This 
edge accumulation is reflected in the Hartree potential as 
the formation of triangular-shaped wells near the edges, 
an effect known in quantum wires prepared by cleaved- 
edge overgrowth 3 ^ as well as from wide two-dimensional 
quantum wells 3 - 3 -. In the non-interacting case, on the 
other hand, the charge density distributions which are 
obtained by using the potential profiles in the left pan- 
els are quite flat, while the electron density oscillations 
between neighboring sites persist. 

For comparison, we also plot the electron density and 
potential of a monolayer GNR of identical width. The 
constant potentials shown in the left panels are classical 
responses of the BGN to the gate voltage correspond- 
ing to the position of the charge neutrality point V c . The 
BGN system can be viewed as two monolayer GNRs con- 
nected in series such that the classical capacitance of the 
BGN is smaller than that of monolayer GNR, resulting 
in larger values of V c . Comparing the averaged electron 
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FIG. 2: Comparison between the single-particle (left parts) 
and the self-consistent Hartree (right parts) calculations for 
different gate voltages. The distance between the gate and 
the BGN is set to d = 100 nm. (a) Induced electron density 
distribution along the y direction (thin lines) and the values 
averaged over six neighbors (bold lines), (b) Potential profiles 
from the one-electron and the Hartree approximation calcula- 
tions, averaged over six neighbors, (c) Band structures. The 
dashed horizontal lines indicate the Fermi level. 



density, one finds that they are very close at the center 
of nanoribbons in both cases, which is consistent with 
the fact that the total induced charge density equals to 
the sum of the densities in each layer, i.e., satisfying the 
relation a = —(fii + p-i) according to electrostatics for bi- 
layer. There are however some small differences close to 
edges which we interpret as a consequence of interlayer 
Coulomb interactions. 

Fig- 121c) shows the dispersion relations of V g = IV 
and 5V for the non-interacting and self-consistent cases. 
Within the single-particle picture, the potential energy 
across the BGN is constant and the dispersion relations 
are shifted upwards almost rigidly as the gate voltages in- 
creases. The energy gaps induced by the gate voltage are 
barely visible. Turning on the electronic interactions, the 
potential profile changes only slightly for V g = IV. How- 
ever, as the gate voltage increases to 5V, the dispersion 
relation gets strongly modified. Anticrossings are gener- 
ated and a significant energy gap is opened between the 
conduction and valence bands, originating not only from 
the potential difference between the layers but also from 
the transverse potential profile. We note in this respect 
that even in monolayer GNRs, the external electric field 
can lead to a gap modulation for semiconductor armchair 
nanoribbons^ and to a gap opening for zigzag ribbons^. 

We continue by considering the average electron den- 
sity and the Hartree potential in each graphenc layer as 
a function of the gate oxide thickness d. As shown in 
Fig. EJa), the electron densities in both layers increase 
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FIG. 3: (a) The average electron density on each graphene 
layer for dielectric thickness d — 1 5nm. Inset: The difference 
of the average density between layers for different dielectric 
thicknesses, (b) The difference of Hartree potentials between 
graphene layers. Solid and dashed lines correspond to the nu- 
merical and analytical calculations from Eq. , respectively. 



approximately linearly V g . In monolayer GNRs, a similar 
linear relation of electron density and gate voltage p oc V g 
has been reported^ This common feature indicates that 
the classical electrostatic contribution predominantly de- 
termines the charge density. Furthermore, the electron 
densities of the two layers of the BGN diverge as the 
gate voltage increases. This electron density difference 
grows monotonously as shown in the inset. The depen- 
dence of the potential differences on d is basically gov- 
erned by the classical component. Fig. [3jb) shows the 
corresponding results of the Hartree potential. The ana- 
lytical expression according to Eq. ^ quantitatively re- 
produces the exact numerical results at the low gate volt- 
ages, especially for larger dielectric thickness d = 50 nm 
and 100 nm. However, deviations emerge for higher gate 
voltages, most pronounced for d = 15 nm. They can 
be understood via the derivation of Eq. which is 

based on the assumption that induced electron density is 
homogeneous and the potential across the ribbon is con- 
stant. In reality, both a larger V g and a smaller d lead to 
the higher electron densities, stronger Coulomb repulsion 
and more pronounced potential inhomogeneities, render- 
ing the classical capacitor approximation less valid. 

We proceed with a detailed study of the effects of 
Coulomb interactions on the band structures, and be- 
gin with looking at the size of the energy gap as a func- 
tion of the gate voltage for different dielectric thicknesses 
as shown in Fig. IHa). In contrast to the 2D bilayer 
graphene, where the band gap increases monotonously 
as the applied electric field is increased^ the energy gap 
in BGNs exhibits a pronounced dip at low gate voltages 
for all values of d. If the single-particle potential is used 
instead, the energy gap increases monotonously with V g 
and follows closely the value of the Hartree potential dif- 
ference shown in Fig. [3Jb), and the dip is absent. This 
implies that the interactions are responsible for the for- 
mation of the dips. To shed light on the origin of this 
suppression of the band gap formation, let us inspect the 
band structures for d = 15nm at several representative 
Fermi energies. For the first two selected values of V g , 
® and (2), the size of gap is proportional to V g . One 
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FIG. 4: (a) The energy gap as a function of the applied gate 
voltage for various dielectric thicknesses d. (b) Representative 
band structures corresponding to the Fermi energies shown in 
(a). The dashed lines indicate the positions of Fermi energy. 



observes that the second energy band of the electrons 
develops an anticrossing with the third one and devel- 
ops a strong curvature around k x — 0. The minimum of 
this band approaches the conduction band edge formed 
by the lowest energy band and takes over the role of 
the band edge at the scenarios depicted in Figs. H^b), 
@ and (3). As a result, the energy gap gets reduced as 
the gate voltage increases and can even approach zero 
for some parameters. We note that a similar effect has 
been discussed by Avctisyan et al., in graphenc tri- and 



quadruple-layer systems based on a different mechanism, 
however J£ As the voltage is increased beyond the dip of 
the energy gap (e.g. V g > 25 for d = 15 nm), the familiar 
Mexican hat shaped dispersion emerges which leads to 
the reduction of the energy gap while the electron-hole 
symmetry vanishes. 

From this evolution of the energy band dispersions, it 
becomes apparent that the Hartree term induces some 
warping of the energy dispersion. We note that in con- 
trast to 2D bilayer systems^ this effect is not due to 
73, since after setting 73 = the band structure changes 
only slightly but the collapse of the band gap persists 
(not shown here). The prerequisite for is the potential 
inhomogeneity generated by the charge accumulation at 
the edges. This is why the collapse of the band gap is ab- 
sent in the 2D cases. Furthermore, the small reduction of 
the energy gap for large V g in the system with d = 15 nm 
is not directly related to warping but is the result of the 
Mexican-hat band formation. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the electron density distribution and 
Hartree potential of a single-gated BGN. Coulomb inter- 
actions are incorporated in a self-consistent way within 
the Hartree approximation. The repelling electrons accu- 
mulate at the BGN edges and induce characteristic dips 
in the transverse potential profile. They change the band 
structures strongly and modify the energy gaps which 
thereby behave qualitatively differently compared to 2D 
bilayer systems, showing an unusual reduction in the sizes 
of the energy gaps at intermediate gate voltages. An an- 
alytical expression for the potential difference between 
the two layers is obtained based on the assumption of 
a classical capacitor and compared with self-consistent 
numerical calculations. They show good agreement at 
small gate voltages. The discrepancies found at higher 
voltages and small dielectric thicknesses are due to the 
assumption of homogeneous charge density and constant 
potential in the analytical model. This becomes invalid 
for high electron densities because the quantum mechan- 
ical effects modify the charge redistribution significantly. 
A disadvantage of the single-gate structure discussed here 
is that due to small capacitance, the Fermi energy is far 
above the neutrality point in the regime of the collapse 
of the band gap, such that this effect is not directly ac- 
cessible experimentally. This shortcoming can concep- 
tually be overcome by performing corresponding simula- 
tions for a double-gate structure, where independent tun- 
ing of the Fermi energy and the interlayer electric field is 
possible 



Acknowledgments 

H.X. and T.H. acknowledge financial support from 
Heinrich-Hcinc-Universitat Diisscldorf and the German 



6 



Academic Exchange Service (DAAD) within the DAAD- (VR) and from the Swedish Foundation for International 
STINT collaborative grant. A.A.S. and I.V.Z. acknowl- Cooperation in Research and Higher Education (STINT) 
edge the support from the Swedish Research Council within the DAAD-STINT collaborative grant. 



Electronic address: hengyi.xu@uni-duesseldorf.de 
' Electronic address: igor.zozoulenko@itn.liu.se 

1 S. V. Morozov, K. S. Novoselov, M. I. Katsnelson, F. 
Schedin, D. C. Elias, J. A. Jaszczak, and A. K. Geim, 
Phys. Rev. Lett. 100, 016602 (2008). 

2 J. Chen, C. Jang, S. Xiao, M. Ishigami, and Michael S. 
Fuhrer, Nature Nanotech. 3, 206 (2008). 

3 M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. 
Phys. 2, 620 (2006). 

4 B. Trauzettel, D. V. Bulaev, D. Loss and G. Burkard, 3, 
192 (2007). 

5 P. San- Jose, E. Prada, E. McCann, and H. Schomerus, 
Phys. Rev. Lett. 102 247204 (2009). 

6 Vitor M. Pereira and A. H. Castro Neto, Phys. Rev. Lett 
103, 046801 (2009). 

7 R. N. Costa Filho, G. A. Farias, and F. M. Peeters, Phys. 
Rev. B 76, 193409 (2007). 

8 E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96 086805 
(2006). 

9 J.B. Oostinga, H.B. Heersche, X. Liu, A.F. Morpurgo, and 
L.M.K. Vandersypen, Nat. Mat., 7, 151 (2007). 

10 Y. Zhang, T. Tang, C. Girir, Z. Hao, M. C. Martin, A. 
Zettl, M. F. Crommie, Y. R. Shen, F. Wang, Nature (Lon- 
don) 459, 820 (2009). 

11 K. F. Mak, C. H. Lui, J. Shan, and T. F. Heinzl, Phys. 
Rev. Lett. 102, 256405 (2009). 

12 F. Xia, D. B. Farmer, Y. Lin, and P. Avouris, Nano. Lett. 
10, 715 (2010). 

13 M. Koshino and E. McCann, Phys. Rev. B 79, 125443 
(2009). 

14 E. McCann, Phys. Rev. B 74, 161403(R) (2006). 

15 A. A. Avetisyan, B. Partoens, and F. M. Peeters, Phys. 
Rev. B 79, 035421 (2009). 

16 A. A. Avetisyan, B. Partoens, and F. M. Peeters, Phys. 
Rev. B 80, 195401 (2009). 

17 H. Min, B. Sahu, S. K. Banerjee, and A. H. MacDonald, 
Phys. Rev. B 75, 155115 (2007). 

18 E. K. Yu, D. A. Stewart, S. Tiwari, Phys Rev. B 77, 195406 
(2008). 

19 P. Cava, M. Lazzeri, A. Marco Saitta, and F. Mauri, Phys. 
Rev. B 79, 165431 (2009). 



20 M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. 
Rev. Lett. 98, 206805 (2007). 

21 X. Wang, Y. Ouyang, X. Li, H. Wang, J. Guo, and H. Dai, 
Phys. Rev. Lett. 100, 206803 (2008). 

22 D. A. Areshkin, D. Gunlycke, and C. T. White, Nano. Lett. 
7, 204 (2007). 

23 K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, Phys. 
Rev. B 59, 8271 (1999). 

24 K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dressel- 
haus, Phys. Rev. B 54, 17954 (1996). 

25 M. Evaldsson, I. V. Zozoulenko, Hengyi Xu, and T. 
Heinzel, Phys. Rev. B 78, 161407(R) (2008). 

26 Hengyi Xu, T. Heinzel, and I. V. Zozoulenko, Phys. Rev. 
B 80, 045308 (2009). 

27 E. R. Mucciolo, A. H. Castro Neto, and C. H. Lewenkopf, 
Phys. Rev. B 79, 075407 (2009). 

28 J. Fernandez-Rossier, J.J. Palacios, and L. Brey, Phys. 
Rev. B 75, 205441 (2007). 

29 N. M. R. Peres, J. N. B. Rodrigues, T. Stauber, and J. M. 
B. Lopes dos Santos, J. Phys. Condens. Matt. 21, 344202 
(2009). 

30 P. G. Silvestrov and K. B. Efetov, Phys. Rev. B 77, 155436 
(2008). 

31 A. A. Shylau, J. W. Klos, and I. V. Zozoulenko, Phys. Rev. 
B 80, 205402 (2009). 

32 S. Ihnatsenka and I. V. Zozoulenko, Phys. Rev. B 74, 



075326 (2006). 



Y. W. Suen, J. Jo, M. B. Santos, L. W. Engel, S. W. Wang, 
and M. Shayegan, Phys. Rev. B 44, 5947 (1991). 

34 A. H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

35 Hengyi Xu, T. Heinzel, M. Evaldsson, and I.V. Zozoulenko, 
Phys. Rev. B 77, 245401 (2008). 

36 C. Ritter, S.S. Makler, and A. Latg, Phys. Rev. B 77, 
195443 (2008). 

37 Y.-W. Son, M.L. Cohen, and S.G. Louie, Nature (London) 
444, 347 (2006). 

38 M. Koshino, New J. Phys. 11, 095010 (2009). 



